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One-loop Feynman Diagrams and Hypergeometric functions. 

Within dimensional regularization the algorithm for analytical evaluation of one-loop multilegs 
Feynman Diagrams has been described a long time ago ||]. It includes two steps: reduction of the 
amplitude to a set of master-integrals [Q] with following analytical evaluation of master-integrals 
via Feynman parameter representation [||]. Due to the appearance of l/£ terms coming from IR 
and/or UV singularities, the NNLO calculation would demand the knowledge of higher terms in 
the £-expansion of the one-loop Feynman Diagrams (see also [Q]). The above mentioned technol- 
ogy is suitable for the evaluation of the finite part of master-integrals [||, However, the direct 
application of this technique for analytical evaluation of the higher-order coefficients in power of 
£ gives rise to complicated results ^, ^ even in a simple kinematic. A perspective approach to 
the construction of analytical coefficients of the £-expansion of one-loop Feynman Diagrams is to 



explore the hypergeometric representation [|10l |11|]. Based on the approach developed in []11|], the 



all-order £-expansion of one-loop self-energy diagrams has been constructed in J = 4 — 2e. More 



examples of hypergeometric representation for one-loop diagrams based on the technique of [10] 



are given in [Il3|,|14j]. However, still by now a systematic way to construct the analytical coefficients 
of the £-expansion for Horn-type hypergeometric functions around rational values of parameters 
does not exist. 

Definition of the liypergeometric function. 

We remind the definition of Horn-type Hypergeometric Functions: 

it is a formal (Laurent) power series in r variables of the following form, 

m,,m2%n,=0 i^b^l ^kbmb + Ok) J 

with Hab,Vah e Q, rj,ak G C and /= {fj,o}. 
The problem under consideration: 

In arbitrary d-dimensional space time, d = 4 — 2e, where £ is the parameter of dimension regu- 
larization the set of discrete parameters /= {7y,5} of hypergeometric function, H{Y;d',z) is 
a linear combination of rational and £-dependent coefficients: = Ao,/t + atS, where Ao,/t and ak 
are arbitrary rational numbers. The Laurent expansion of the hypergeometric function around the 
integer value <i = 4, is called "construction of the analytical coefficients of £-expansion" of the 
function: 

H{Ao + ae-z) = H{Ao;z) + £ £^/J;(z) , (2) 



where symbolically. 



hj{^) = ^H{A;z) 



(3) 

A=Ao 



The goal is to write the coefficient functions h; in terms of known functions, suitable for numer- 



ical evaluation [15], or to describe all analytical properties of hj, treating them as a new class of 
functions. 

Existing algorithms: 

The first systematic algorithms for the construction of higher order coefficients of the £-expansion 
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of multivariable hypergeometric functions around integer values of parameters were suggested in 
[16]. In Ref. [|l^], the special set of rational values of parameters, the so called "zero-balance case" 
was analyzed. However, the partial results of [17] beyond the zero-balance case are in contradiction 
with partial results of Ref. [18]. 
Our method: 

In a series of papers [19, 20, ^ it was shown that for the hypergeometric function of one variable, 
pFp-i, the analytical coefficients of the £-expansion can be constructed via an explicit solution of 
differential equations for coefficients functions hj{z). Using Eq. (|3|) it is easy to show that the 
coefficients hj{z) satisfy the following linear system of (partial) differential equations (PDE): 



£p._//(A;?)=0: 



d_ 



lPz§//(A;z)=0 

. L 



A=Ao 



A=Ao 



A=Ao 



(4) 



A=Ao 



When are non-homogeneous PDE solvable in terms of multiple polylogarithms? 

We are interested in the question under which conditions the functions hj{z), solutions of Eq. (Q), 
are expressible in terms of multiple polylogarithms [2^, 23, 25], or generalized iterated inte- 
grals, defined as: 



Giz;Rk,Rk u- -,Ri 



dt 



■I{t;Rk-u---,Ri) 



dtk /■'* dtk- 



Rk{t) Jo Rk-\{t) Jo R\{ti 



^2 dt^ 



(5) 



/o Rk{t 

where Rk{t) are some rational functions. Multiple polylogarithms correspond to Rk{t) = t — a^- 
When is the system of PDE with non-zero non-homogeneous part solvable in terms of (generalized) 
multiple polylogarithms? Our algorithm includes the following steps: 

• Factorization: the differential operator(s) after £-expansion are factorisable into product of 
differential operators of the first order; 

• Linear parametrization to all orders in e; 

• The non-homogeneous part belongs to the class of functions of the special type (see below). 
Example I 

Let us consider the differential equation related to the hypergeometric function pFp^\ [20, 21]: 



l^Pk{z;e) 

k=0 



H{z;e)=F{z;e) , 

where Pk{z;s) and F{z;£) are rational functions or iterated integrals over a rational 1-form: 



Pkiz;£) 



Ujiz-aj-pje) 



F(z;e) 



dt 



Ujit- 



vje) 



t-aU,-[t-Mr-Nre) 



Uriz-Ar-BrS) 

We are looking for a solution of Eq. (^ of the following form: H{z;e) = T.'j=ohj{z)s^ . 
Factorization. Factorization of differential operators after £-expansion means the following: 

k 



(6) 



(7) 



k=0 



in; 

r=0 



lr<P 

k=l 



RkAz)-r + QkAz) 

dz 
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where Rk,r{z) and Qk,r{z) are some rational functions. 

Linear parametrization. Let us consider as illustration the following differential equation 



Ri{z)- + Qiiz) 
dz 



R2{Z)- + Q2{Z) 

dz 



h{z)=F{z) 



(8) 



Its iterated solution is: 



Kz) 



dh 



R2{h) 



exp 



-Jo W"^^ 



'3 dt 



Ry{h) 



exp 



(9) 



In accordance with Eq. this iterated integral can be written as multiple polylogarithm, if there 



is a new variable ^ : ^ = ^{t), converting this expression into ratio of polynomials [20|: 

dt 



(10) 



where Mi,Ni,Ki,Li are polynomial functions. The existence of such a parametrization we called 
linear parametrization. 

When does such parametrization exist? To answer to this question, the non-homogeneous part 
of differential equation Eq. (^ should satisfy the system of linear PDE with Factorization and 
Linear parametrization in each order of e: 



(11) 



where /', ,5';, T are rational functions. 
Multivariable generalization 

Generalization of this technique for the Horn-type hypergeometric functions is straightforward: 
1. Convert the system of linear PDE with polynomial coefficients into Pfaff form: 

£P/.^(a;z)^i^(a;z) = ^ |cfytW;(?) = ^\j{z)(Oj{z)dzk , dr [dk(Oi{z)] = o| . 



2. Find the values of parameters when the last system of linear PDE can be converted into 
triangular form and when Factorization is vahd. 



3. Find a linear parametrization: validity of Eq.( 10) for each variable 



Simplification of the procedure of Factorization. 

To simplify the procedure of Factorization of differential operators, the following trick is very 
useful. Any Horn type hypergeometric function, defined by Eq. (|l|), satisfies the system of linear 
PDE with polynomial coefficients: 



Pliz)^H{J;z)=0 



(12) 
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where ^ = and -P^(z) are polynomial. Moreover, there are linear differential operators 



that change the value of each parameter by ± 1 : 



(13) 



In accordance with p6[], the differential operators inverse to the operators defined by Eq. ( |13[ ) can 
be constructed. Applying direct/inverse differential operators to the hypergeometric function the 
values of parameters can be changed by an arbitrary integer numbers: 



Qo{z)HiJ + m;z) = £ Qjiz)dyH{J;z) 

7=0 



(14) 



where m is a set of integers and Qq{z) and Qj are polynomials and r is the holonomic rank of 
system (12). More details are given in [p7|]. 
Example II: Ft, hypergeometric function 

In |2l]] we apphed our algorithm for the construction of the e-expansion of Fi and F3 Appell 



hypergeometric functions. The results of the e expansion for Fi— [28] and Ft,— functions [21] 
around integer values of parameters are in agreement with results of [^]. Let us consider the 
£-expansion of the Appell hypergeometric function Ft, around rational values of parameters. The 
£-expansion around this set of parameters does not follow from the algorithms described in [|l6[ p^ ]. 
Let us consider the Appell hypergeometric function F3 : 



Ft, hai£, \-a2e, — h^i£,— +^2£,1 \-ce;x,y 

q q q q q 



00 CO 

II 

m=On=0 



m! n\ 



(15) 



m+n 



where {pi,rj,p,q} are integers. Applying our technology step-by-step, we find that the system 
of linear PDE for the coefficient functions is factorisable and has a triangular form only when 
pjrj = for J = 1,2. After that, the original system of linear PDE with polynomial coefficients is 
transformed into a linear system of PDE with algebraic coefficients. To convert this system into a 
class of linear PDE, the linear parametrization should exist simultaneously for the each element of 
the singular locus of Ft, : 



{x}U{l-x}U{3;}U{l-3'}U{x3;-x-3;} 
as well as for the auxiliary functions Hj, j = 1,2,3 defined as 



Hi{x) = (-1) 



SI+S2 

H3{x,y) = (-1) — 



(x- lY'+p 



//2(x) = (-l)t 



(3'-l)^= 



+p 



{xy - X - yY^+'^+P 

where sj = pj + rj and j = 1,2. We find that the linear parametrization exists when: 



(16) 



(17) 
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• The functions Hj are constant polynomial: si = S2 = p = 0. It corresponds to [16|. 

• One of thi^ee functions Hj are equal to 1: i'l = 0,^2 = —p, and (1 -f-)- 2). 

• Two of three functions Hj coincide: si ^ 0,S2 = p = 0, and (1 -f-)- 2). 



Unfortunately, for another physically interesting set of parameters [14|, we failed to rewrite the 



iterated integral in terms of multiple polylogarithms. For example, for = ^2 = and p 0, 
the statement about the existence of a linear parametrization is equivalent to the existence of three 
different (rational) polynomial functions of two variables P\{x,y),P2{x,y) and P3(x,j), such that 



(18) 



where q is integer and q > 2. To our knowledge, this equation has a solution only in the class 
of elliptic functions. However, the finite part of the Fs-function with this set of parameters is 



expressible in terms of polylogarithms []30|]. 

As result we got, that only for the two cases: 



^3 ( h + — +ai£,l2 + a2£,l3 + bi£,l4 + b2£,l5 + — +ce;x,y 
q q 

F3 [ h + ^ +ai£,l2 + a2£,h + bi£,l4 + b2£,l5+c£;x,y 



(19) 
(20) 



where Ij,pi,q are integers, the analytical coefficients of the e-expansion of F3 hypergeometric 



function are explicitly expressible in terms of multiple polylogarithms []24|]. 
Hypergeometric Functions vs. Feynman Diagram 

These two building blocks. Factorization and Linear parametrization, are sufficient to rewrite an 
iterative solution of system of linear PDE in terms of multiple polylogarithms. It is a cornerstone 
of all modem multiloop analytical evaluations of master-integrals in QCD and our results are in 



full agreement with available QCD calculations [^, [33[ ]. 
Conclusion: 

The algorithm described in [|^, ^ 21 ] has been applied to the construction of the analytical co- 



efficients of £ -expansion of Horn-type hypergeometric functions of two variables [ ]34| ] as well as 



Mellin-Bames integrals [35|. In particular, we analyzed the following linear system: 



Uodn(o(z;£) = { Uidi2+Piei+P2d2+Po U(z;£) , 



Tod22(0{z;£) = { Tidi2+Rldi+R2d2+R3 >«(?;£) , 



(21) 



where z = (zi,Z2) are independent variables, dj = zjd^. , j = 1,2 , and = 0,- • • • The func- 
tions Go = {Uq, Tq,Ui,Ti} are polynomial of variables zi and Z2- 



(22) 
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all other function Eo = {Pi,Rj} are polynomial of variables zi,Z2 and £: 

^0= £ YijAzie' , (23) 

ij,k=0 

and Gij and jij-k ai'e rational. The analytical structure of the coefficients h^p (z) defined via Laurent 
expansion around e = of the solution (0^''\z',£) of this system have been analyzed 



«W(z;£) = £/iW(z)£^ (24) 



for the physically interesting set of parameters [14|. In particular, the hypergeometric functions 



considered in ^ correspond to a system of linear PDE ( pl| ) with polynomial coefficients and 
with singularity locus 

L := {zi } U {Uo} U {zi} U {ro} U {UoTo -UyTi} , (25) 

where 

Ui{z;a) = aoj + aijZi+a2jZ2 , Tj{z;a) = bijzi +^2jZ2 + ^3.;ZiZ2 , (26) 

ij = 1,2 and akj,hj G {0,±1}. See also 

The £-expansion around rational values of parameters with one unbalanced rational parameter, 
corresponds to a system of linear PDE with algebraic or elliptic coefficients. Imposing only Fac- 
torization conditions gives rise to iterative integrals with algebraic functions, that in general, are 
not expressible in terms of multiple polylogarithms. Only when additional Linear parametrization 
conditions are valid, we are able to rewrite the results of the integration in terms of 2-dimensional 
polylogarithms [|^. The Linear parametrization should exist simultaneously for the each element 
of singular locus, Eq. (25), of the differential system Eq. (|1|) and for algebraic functions defined as 



q-roots of ratios of elements of L: (L,/(l — L;))« and/or ((L,Ly)/(L,+Lj — L,Ly))« (see Eq. (17)). 
It is in agreement with the one- variable case analysed in [pO|]. 

We got, see also [pO|], that even when the finite part of a hypergeometric function is expressible 
in terms of multiple polylogarithms (existence of Liouvillian solution of a linear system of PDE 
in J = 4) it does not follow that higher order terms of the £-expansion are expressible in terms of 
multiple polylogarithms, too. 
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